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A theoretical model for understanding photomagnetism in the heptanuclear complex 
[Mo(IV)(CN)2(CN — CuL)e] 8+ perchlorate is developed. It is a many-body model involving the 
active orbitals on the transition metal ions. The model is exactly solved using a valence bond 
approach. The ground state solution of the model is highly degenerate and is spanned by five 
S=0 states, nine S=l states, five S=2 states and one S=3 state. The orbital occupancies in all 
these states correspond to six Cu(II) ions and one diamagnetic Mo(IV) ion. The optically excited 
charge-transfer (CT) state in each spin sector occur at nearly the same excitation energy of 2.993 
eV for the physically reasonable parameter values. The degeneracy of the CT states is largest in 
the S=3 sector and so is the transition dipole moment from the ground state to these excited states. 
Thus laser irradiation with light of this energy results in most intense absorption in the S=3 sector. 
The life-time of the S=3 excited states is also expected to be the largest as the number of states 
below that energy is very sparse in this spin sector when compared to other spin sectors. These twin 
features of our model explain the observed photomagnetism in the [Mo(IV)(C N)2(C N — CuL)e] 8+ 
complex. 
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I. INTRODUCTION 

In the expanding field of molecular magnetism, pho- 
tomagnetism has considerable import as it affords the 
possibility of magnetic switching with the aid of light. In 
the photomagnetic systems that we have studied here, 
the system before shining light is essentially paramag- 
netic, with negligible exchange interaction between the 
magnetic ions. Shining a burst of light of appropriate 
wavelength takes the system to a state in which there is 
a superexchange interaction between the paramagnetic 
ions. The photo-induced high-spin state is long lived at 
low temperatures, however, warming the system takes it 
back to the paramagnetic ground state [l|. There also 
exist extended systems in which the photoinduced ex- 
change mechanism gives rise to magnetic ordering; this 
unusual photo-effect is called the Photo Induced Phase 
Transition (PIPT) 0. 

The above phenomenon is quite distinct from the spin 
state transition observed in many inorganic complexes, 
wherein shining light leads to change in the ratio of 
the population of the low-spin (LS) state to that in 
the high-spin (HS) state. This phenomenon called the 
Light Induced Excited Spin State Trapping (LIESST), 
has been extensively studied in spin crossover complexes 
. These systems lie very close to the spin crossover 
point in the Tanabe-Sugano diagram 4]. When the sys- 
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tern is excited by light, the non-radiative decay from the 
excited state occurs through intersystem crossings to the 
two nearly degenerate low energy states of different spin 
5]. To recover the pure low-spin ground state, it would 
become necessary to warm the system and cool it to low 
temperatures. This phenomenon does not involve any 
magnetic exchange interactions between the transition 
metal ions. 

Photomagnetism was experimentally first 
demonstrated in Fe-Co prussian blue complex 
K .2Coi.4Fe(CN) 6 -6.9H 2 O. This system is ferrimag- 
netic with a T c of 16K Upon irradiation, the T c 
increased to 19K. Besides, at fixed magnetic field and 
fixed temperature, the net magnetization increased upon 
irradiation. This increase in T c and magnetization are 
attributed to electron transfer from the diamagnetic 
Fe(II) to the diamagnetic Co(III) brought about by 
irradiation resulting in Fell 1 1) ion with spin 1/2 and 
Co(II) ion with spin 3/2 These paramagnetic ions, 
besides contributing to an increase magnetization, also 
establish pathways for exchange which were blocked in 
the diamagnetic state, thereby increasing the T c . 

The system we have studied theoretically is the Mo — 
Cu molecular system. While photoinduced ferromag- 
netism has been observed in extended Mo — Cu systems 
, we focus our attention in this paper on the molecular 
system [Mo(IV)(CN) 2 (CN - Cu(II)L) 6 } 8+ with per- 
chlorate as the counter ions Q. In the ground state, 
each molecule consists of six noninteracting spins, de- 
rived from the six Cu(II) ions; the Mo(IV) ion is dia- 
magnetic. The magnetic behaviour of this molecule in 
the ground state corresponds to that of a system with 
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FIG. 2: Topology of the MoCuq system with Mo ion at the 
center and Cu ions at the vertices of a regular hexagon. All 
the Mo — Cu distances are 5.25 Afrom Mo. 



II. MICROSCOPIC MODEL 



FIG. 1: Temperature dependence of XmT, (o) before irradi- 
ation, (•) after irradiation. 

six independent spin-1/2 objects per formula unit with 
XmT value of 2.48 cm 3 K mol -1 at 20K. On irradiation 
with blue light (406-415 nm) at 10K, the xmT value in- 
creases to 4.8 cm 3 K mol -1 . On turning off the light and 
warming the system, the xmT value peaks at 5 cm 3 K 
mol -1 and decreases thereafter gradually to a room tem- 
perature value of 2.6 cm 3 K mol -1 . The experimental 
data after irradiation is in agreement with 75% of the 
molecules in S=3 state with the remaining in the state 
corresponding to six isolated S=l/2 spins per molecule. 
Fig. ^shows the XmT vs T plot of the system before and 
after irradiation. 

In this paper, we develop a microscopic quantum 
many-body model to understand the photomagnetic be- 
haviour of this system. In the next section, we describe 
our model and the computational scheme in detail. In 
the third section we present the results from our model. 
We show that the eigenvalue spectrum of the model can 
be classified according to the occupancy of the various or- 
bitals. The low-energy eigenvalue spectrum corresponds 
to one doubly occupied and one empty Mo orbital and 
singly occupied Cu orbitals with all the spin states in 
this sector being degenerate. In each of these spin sec- 
tors, we have several states to which approximately 3eV 
optical excitation has nonzero transition dipole moment. 
For a reasonable set of parameters the transition dipoles 
to the S=3 excited states are the largest. We also find 
that the S=3 density of states (DoS) is an order of mag- 
nitude smaller than the DoS in other spin sectors. This 
would lead to long life-time for the S=3 optical excita- 
tions resulting in the system behaving like a S=3 molec- 
ular complex after irradiation. 



The topology of the arrangement of the transition 
metal ions in the complex is hexagonal with the Mo ion 
at the center of the hexagon and the Cu ions at the ver- 
tices (Fig. |2l . The molecular geometry is such that there 
is no electron transfer pathway which directly connects 
the copper ions, hence all exchange interactions between 
the Cu ions occur through the central Mo ion. Low- 
symmetry of the complex around the Mo ion implies that 
the degeneracy of the 4<i orbitals on the Mo ion is lifted. 
This allows one to model photomagnctism in this sys- 
tem by considering only two low-lying Ad orbitals on Mo 
and the partly occupied 3c? orbital on each Cu. If we 
restrict the number of Mo orbitals to just one, no sig- 
nificant exchange interaction between the paramagnetic 
sites is possible since intra-atomic direct exchange inter- 
action would be absent. Each Cu(II) ion contributes one 
electron to the system while the Mo(IV) ion contributes 
two electrons. Thus, we would be dealing with an open 
shell system of eight electrons in eight orbitals. 



A. Model parameters 

The essential parameters of the model are the effec- 
tive transfer parameter (the hopping integral) between 
the Cu orbital and each of the two orbitals on Mo. It is 
assumed that the hopping integral, t, between a Cu or- 
bital and any of the two active Mo orbitals is the same. 
The parameter A defines the orbital energy of the lower 
energy orbital on the Mo site, with reference to the or- 
bital energy of the Cu d-orbital taken to be zero. The 
energy splitting between the two orbitals on the Mo site 
is given by the parameter 5. The quantities t, A and 
5 are all parameters of the non-interacting part of the 
Hamiltonian of the system. 

The two-electron integrals, which are important for 
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magnetism and should appear in the interacting part of 
the model Hamiltonian are the intra-orbital, the on-site 
inter-orbital and inter-site inter-orbital electron repulsion 
integrals. The first of these is given by, 

U i = [ii\i{\ = J f^*(l) ( f> i a)^*(2)4> i {2)d 3 r 1 d 3 r 2 

(1) 

where U l is the familiar Hubbard parameter. We have as- 
sumed that the two Mo orbitals have the same Hubbard 
parameter. In cases where a given site has more than one 
orbital, we need to include inter-orbital Coulomb and ex- 
change electron repulsion integrals in the model. There 
are three kinds of inter-orbital repulsion integrals, with 
two orbitals. These are denoted as U m , X l% and W" 
and are given by, 

U u ' = = f f 0*(l)0. i (l)-|-0:,(2)^(2)d 3 r 1 d 3 r 2 

<j>* (1)4' (1)— 0*(2)</v (2)d 3 nd 3 r 2 

W u ' = [ii\ii'} = J /"#(l)&(l)i#(2)&, (2)^^(2) 

where we have used the charge cloud notation of 
chemists' to define the two electron integrals both here 
and elsewhere in the paper HJ. The single occupancy of 
the orbitals 4>i or 4>i' , is favoured over double occupancy 
of either, when the inter-orbital repulsion U li is smaller 
than the Hubbard parameter U l or U l for either of the 
orbitals <pi or <pi> . The X integral is the exchange integral 
which favours parallel alignment of spins in the orbitals 
4>i and 4>i> when the orbitals are singly occupied, lead- 
ing to the familiar Hund's rule. The W zt term does not 
play a significant role in magnetism but is included in the 
model for the sake of consistency since it is larger than 
the exchange integral. In our model, these terms arise 
only in the case of the Mo atom as it contributes two 
active orbitals to the model, say 1 and 2, thus we have 
U 12 = [11 1 22], W = [11|12] and X = [12|12] We 
also include the extended range electron correlation cor- 
responding to repulsion between charge densities located 
on different sites, within the zero differential overlap ap- 
proximation, via Ohno interpolation |ll| . The model 
Hamiltonian can now be written as 
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Here, H.c. stands for Hermitian conjugate; our number- 
ing scheme for the orbitals is such that orbitals 1 and 2 
are on Mo site and orbitals 3 through 8 are each located 
on the copper sites and 



Kj = 4<A>! = a\ *a ha (4) 

cr cr 

The inter-site interaction V%j is parametrized as 



Vij = 14.397 



28.794 

w + w 



(•5) 



—2{Ei2 + E 2 i] 



with the distance is in A and energies are in eV. The 
geometry of the cluster MoCuq is taken to be a regular 
hexagon with a Mo — Cu distance of 5.25A. The topology 
of the complex used in modeling the system is shown 
schematically in Fig. 



B. Computational details 

The model Hamiltonian spans a finite dimensional 
Hilbert space corresponding to eight electrons in eight 
orbitals. Since the Hamiltonian conserves total spin, the 
Hilbert space can be further subdivided into spaces of 
total spin 0, 1, 2, 3 and 4 whose dimensions are 1764, 
2352, 720, 63 and 1 respectively. The basis states with 
a given total spin can be generated using a valence bond 
approach of explicit spin pairings |l2j . For example, if 
singly occupied sites i and j are spin paired, then a line 
is drawn between sites i and j in the VB diagram to in- 
dicate the spin pairing, («i/3j — /3iOj)/v2, i < j. We say 
that the line begins at site i and ends at site j. If the 
spins at sites ii,i2 • ■ ■ ii are not paired, we pass an arrow 
through these sites in the VB diagram and this denotes 
the state a^a^ ■■■ai l . Because the Hamiltonian con- 
serves Sf otal besides S* 2 , it is sufficient to work in the 
subspace Ms = S. The complete and linearly indepen- 
dent set of VB states in each spin space can be obtained 
by taking recourse to modified Rumer-Pauling rules . 
Some typical VB diagrams are shown in Fig. 01 Each 
orbital in the VB picture has one of four possibilities; 
the orbital is (i) empty, (ii) doubly occupied, (iii) a line 
begins at the orbital or (iv) a line ends at the orbital. 
It is possible to associate these four possibilities of an 
orbital in a VB diagram with the four states of two bits 
and stored as sixteen bit integers in an ascending order. 
In any given spin space, the effect of the term £!y in the 
Hamiltonian on a basis state is to alter the orbital occu- 
pancy of orbitals i and j (subject to Pauli principle) and 
pair the spins in the orbitals that were involved to yield 
a new VB digram with a fixed amplitude. If the new 
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lence bond basis to the constant Ms basis and carry out 
all these calculations. 



III. RESULTS AND DISCUSSION 
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FIG. 3: Representative VB diagrams in (a) S=0, (b) S=l, 
(c) S=2 and (d) S~3 subspaces. Dots (•) represent empty 
sites and crosses (x) represent doubly occupied sites. In each 
VB diagram, numbers 1 through 8 are the site indices and the 
values at the bottom against each site gives the bit state of 
the bits representing the site in the integer I. The bit-pattern 
and the value of the integer I, representing the VB diagram 
are also shown. 



VB diagram violates Rumor- Pauling rules, it is trivially 
possible to express it as a linear combination of the basis 
VB states. 

Using the above procedure, the Hamiltonian matrix 
can be set-up in the chosen total spin sector. The result- 
ing Hamiltonian matrix is sparse (since the number of 
terms in the Hamiltonian is far smaller than the dimen- 
sion of the VB complete space) and nonsymmetric matrix 
and can be partially diagonalized to obtain a few of the 
low-lying eigenstate s, using Rettrups modification of the 
Davidson algorithm (14lll5l |. In some cases, full diagonal- 
ization of the matrices is resorted to obtain the complete 
eigenvalue spectrum. In our case, since the dimensional- 
ities of the different subspaces are fairly small, we have 
resorted to complete diagonalization of the Hamiltonian 
matrices in each spin sector. 

In the present problem, the quantities of interest are 
the optical gaps in each spin sector and the transition 
dipoles for these transitions. To understand the nature 
of the excited states, it is also important to know the spin 
and charge density distributions in these states. Comput- 
ing many of these quantities is simpler in the constant 
Ms basis. Hence, we transform eigenstates in the va- 



We have extensively explored the properties of the 
model over a wide range of parameters, specifically, we 
have studied the model at discrete values of the parame- 
ters in the following ranges: 4 eV < U Cu < 6 eV, 0.5 eV 

< U Cu - U Mo < 2.0 eV, -0.2 eV < A < -0.5 eV, 0.0 eV 

< S < 0.3 eV, 0.5 eV < u Mo -U M °- Mo < 1.5 eV, 0.5 eV 

< W < 1.5 eV with X < W/2. However, we report our 
results for the model parameters values of U Cu =6.0 eV, 
C/ Mo =4.0 eV, A=-0.4 eV, <5=0.0 eV, U Mo ~ Mo = 3.5 eV, 
IT=0.75 eV and X—0.3 eV since it gives a reasonable 
value for the optical gap of the system. Our assump- 
tion of U Cu being higher than JJ Mo is justified since 3d 
orbitals are more compact than the Ad orbitals. The in- 
terelectron repulsion integral involving the d orbitals on 
Mo, (£/ M °- Mo ) i s slightly less than the intra-orbital re- 
pulsions. We expect the integral W to be larger than X, 
and both are expected to be significantly smaller than 
the U type repulsion integrals that involve repulsion be- 
tween two electron densities, each one located in a single 
orbital. The orbital energy of the Mo Ad orbitals could 
be slightly smaller than the Cu 3d orbital energy since 
Mo has a higher atomic number. Therefore, we have 
assumed A to be small and negative. These values of 
the interaction parameter are reasonable as found from 
theoretical studies on transition metal oxides in the con- 
text of high-T c superconductors .16]. For these model 
parameters values, transfer integrals in the range 6xl0 -4 
cV < t < 0.01 eV fit the experimental x T vs T data. This 
is because, the various spin states with all Copper ions in 
+2 oxidation state and Molybdenum ion in +4 oxidation 
state are degenerate to within about t 2 /U, where U is 
an average Hubbard parameter (~ 5 eV). The near de- 
generacy of these " covalent" spin states ensures that the 
behaviour for T > t 2 /U corresponds to that of six free 
spins on the Cu(II) ions. We have fixed the value of t 
at 0.01 eV to ensure that there is a significant transition 
dipole between the low-lying spin state and its optically 
coupled state. 



A. Low-energy spectrum of the model 

To understand the eigenvalue spectrum of the model, it 
is instructive to analyze the t — behaviour. In this case, 
the eigenstates of the model can be completely classified 
based on the site occupancies. The only off-diagonal con- 
tribution that mixes the states with different occupancies 
arises from the W and the X terms involving the Mo or- 
bitals. These terms, however, conserve the total number 
of electrons on Mo. When we have the Ad 2 configuration 
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FIG. 4: Theoretical fit of temperature dependence of XmT for (a) t = 0.01 eV and (b) t = 0.00065 eV. Error in the fit R 
calculated from R = J2 l (Xe XP t(T l ) - Xcaic{T % )f /E 4 Xex P t(T<) 2 for (a) 1 x HT 5 and (b) 1.35 x 10" 5 . 



on Mo and all the copper ions in 3d 9 configuration, wc 
have a total of 384 states. This corresponds to six ways 
of occupying two electrons in the two 4d orbitals and 
sixty four spin configurations of the six spins on the six 
Cu(II) ions. For the model employed by us, the lowest 
energy state corresponds to the lower of the orbitals on 
Mo being doubly occupied, for U Mo < (U M °- Mo + S), 
ignoring the off-diagonal contributions from the W and 
X terms. The effect of the off-diagonal W and X terms is 
(i) to favour double occupancy on the lower energy Mo 
orbital and (ii) to favour parallel spin alignment when 
these orbitals are singly occupied. Corresponding to this 
occupancy, five states with total spin S=0, nine states 
with total spin S=l, five states with total S=2 and one 
state with total spin S=3 are all degenerate. There is 
only one S=4 state which is 0.671 eV above the ground 
state and corresponds to single occupancy of each orbital 
with all spins aligned in parallel. These states span all 
the 64 spin configurations corresponding to six spin-1/2 
objects; since all these states are degenerate to within a 
few degree K, the paramagnetic behaviour observed in 
the system before irradiation is reproduced for small val- 
ues of t. The theoretical fit to experimental data for two 
different values of the transfer integral is shown in Fig. 

El 

The next set of states corresponds to single occupancy 
of each of the Mo orbitals. If the spins of the electrons 
in these orbitals are parallel, the states will have a lower 
energy than if they are antiparallel. There are in all 
256 such states of which 128 states have parallel spin 
alignment and another 128 states have antiparallel spin 
alignment. The remaining sixty four of the 384 states 
have a double occupancy in the higher energy Mo orbital. 



These 384 states form the low energy spectrum when U 
is significantly larger than (U Mo — A) so that the intersite 
interactions do not alter the picture. The most important 
point to note is that since the low-energy states all have 
identical charge distribution, the transition dipole matrix 
elements between these states are zero and this result 
holds even when t ^ 0. 



B. Charge- Transfer Excited States 

The higher energy excited states of the system in- 
volve charge-transfer (CT) from the Mo(IV) ion to 
Cu(II) resulting in a state with a large weight for the 
Mo(V) — Cu(I) configurations. These states are found 
at about 3 eV above the ground state. These excita- 
tions have nonzero transition dipole to the ground state. 
One of the most important features of this excitation is 
that, in the excited state there is considerable weight 
for the basis state with singly occupied Mo site. The 
singly occupied Mo(V) site favours virtual electron trans- 
fer from the Cu(II) sites, of electrons which have the 
same spin orientation as the unpaired MoiV) electron, 
because of the Hund's rule. The virtual transfer of an 
electron from Cu(II) with its spin orientation same as 
that already present on Mo(V) is favoured as long as (i) 
the (i-orbitals on Mo are nearly degenerate and (ii) the 
inter-orbital d — d repulsion, on Mo is weaker than the 
intra-orbital repulsion (jj M °- Mo < U Mo ). This feature 
contributes to stabilizing the high-spin CT excitation rel- 
ative to the low-spin excitations. In Fig we show the 
possible pathways for superexchange in this state. The 
factor which favours stabilization of the low-spin state 
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FIG. 5: The virtual excited state resulting from the electron 
transfer shown in (a) has a lower energy than those shown in 
(b) and (c), because of the Hund's rule. The superexchange 
process in (a) favours S=3 state, and in (b) an S=2 state. 
The process (b) is preferred over (a) if the splitting of the 
orbitals on the Mo site as well as jj m °- Mo a re large. In all 
the figures, the two orbitals in the center are the Mo orbitals 
and the six peripheral orbitals are the Cu orbitals. 



is the large number of pathways that exist for the delo- 
calization of electrons in the low-spin state, compared to 
the high-spin states. In our case, there are 1764 singlets 
and 2352 triplets while there are only 63 S=3 states. The 
splitting of the states with different spins in the CT state 
is thus a competition between direct exchange which sta- 
bilizes the high-spin state and the kinetic exchange which 
favours the low-spin state due to the availability of a 
larger phase space for electron delocalization. This im- 
plies that the lowest-energy CT excited state could have 
any spin, S, in the range < S < 3. 

If the splitting of the Mod orbitals is large or U Mo ~ Mo 
is large, then the intermediate virtual state with a dou- 
bly occupied d orbital on Mo is favoured, resulting in 
an antiferromagnetic interaction between the Cu(II) ion 
and the Mo spin (FigjSJ. This factor would contribute 
to a lowering of the S=2 CT excited state relative to CT 
excited states with spin S^2. However, the final outcome 
depends on both the direct exchange and the kinetic ex- 




FIG. 6: Density of states versus energy in S=0,l,2 and 3 
subspaces. There is only one state in S=4 sector at 0.6709 eV 
from the ground state (not shown in the figure). 



change contributions. Thus, it is again difficult to predict 
which of the spin states is a low-energy CT excited state. 
We see all the above scenarios in our calculations, when 
we widely vary the model parameters. 

The above picture is slightly altered by the long-range 
interactions in the model. These interactions depend 
upon the occupancies of the Mo site as well as all the 
Cu sites and are significant since the molecular complex 
is insulating and screening of these interactions is not 
significant. Together with a nonzero transfer integral, 
the long-range interactions complicate the picture while 
retaining the general physical features outlined above. 

For the model parameter for which we are reporting 
the results, we show the histogram giving the number 
of states as a function of energy in Fig. HJ for the full 
spectrum of the model in each total spin sector. We find 
CT excitations at about 2.993 eV in all the spin sectors. 
The energies and magnitude of the transition dipole vec- 
tors from the ground state are given in Table |U We also 
present in Table [I] the degeneracy of states as well as the 
charge and spin density on the Mo site. We note that 
in all the states located at this energy independent of 
the total spin, the charge density on the Mo site is one. 
What is interesting is that while the energies of these 
states as well as the charge density on the Mo site are all 
the same, the transition dipole moment is largest for the 
S=3 states. This fact, coupled with the high degeneracy 
of the state implies that the S—3 states have the largest 
absorption to the CT excited state. 
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TABLE I: Energy gaps (in eV) to the charge-transfer excited states from the lowest state in different spin sectors and the 
magnitude of the corresponding transition dipole moments in 1CT 3 Debye. Also presented in the table are the total spin and 
charge densities on the Molybdenum site in the excited state. 



Spin sector 


Energy gap (eV) 


Degeneracy 


Transition dipole 
(xl(T 3 Debye) 


< S?° > 


< UMo > 


S=0 


2.9921 


6 


7.0 


0.0000 


1.0000 




2.9923 


4 


4.1 


0.0000 


1.0000 




2.9920 


4 


6.2 


-0.2498 


1.0000 


S=l 


2.9923 


8 


5.3 


0.4986 


1.0000 




2.9924 


2 


3.4 


0.5000 


1.0000 




2.9921 


2 


8.3 


-0.1596 


1.0000 


S=2 


2.9923 


2 


12.1 


0.5000 


1.0000 




2.9924 


4 


7.8 


0.5000 


1.0000 


S=3 


2.9925 


5 


16.0 


0.5000 


1.0000 



An important difference between the S=3 eigenvalue 
spectrum and eigenvalue spectra of other spin states is 
that the DoS in the former is very dispersed and the 
states are far apart in energy, when not degenerate. This 
has implications for the life-time of the excited states. 
In the S=0, 1 and 2 spaces, there are an order of mag- 
nitude more states near the CT excitation gap of 2.993 
eV compared with the S=3 states at this energy. Thus, 
other spin states have pathways for internal conversion 
leading to rapid de-excitation of the CT state. This is 
less likely in the S=3 state. For this reason, the S—3 
excited state has a longer life-time and can be observed 
in the photomagnetism experiments. One feature that is 
however missed in our studies is the computation of the 
equilibrium geometry in the CT excited state specifically 
around the Cu{I) site that may play a crucial role in in- 
creasing the life-time of the state. This is rather hard to 
incorporate in our model wherein the ligand atoms are 
completely neglected. It is possible that the CT excited 
state with S=3 is more distorted (relative to the ground 
state) than similar states with different total spin. This 
could further enhance the life-time of the S=3 CT exci- 
tation. 

For the sake of completeness, it is worth mentioning 
that there are also different kinds of charge transfer ex- 
cited states at much higher energies corresponding to 
electron transfer from the Cu(II) site to the Mo(IV) 
site, charge transfer between Cu(II) sites and so forth for 
which the transition dipole moments are nonzero. How- 
ever, these states occur at a much higher energy and are 
not excited by laser light at 3eV. 

The model developed in this paper is specific to the 
isolated molecular complex we have studied. However, 
this model also throws light on the mechanism leading 
to the observed ferromagnetism in the PIPT systems. 
While the degeneracy of the different optically excited 
spin states is not significantly affected by the exchange 
mechanism operative in the optically excited state, the 
long life-time of the high-spin state could be responsible 



for the observed ferromagnetism in extended systems. 



IV. CONCLUSIONS 

In this paper, we have introduced a model for 
the photomagnetism in the hcptanuclear complex, 
[Mo(IV)(CN) 2 {CN - CuL) e ] s + perchloratc which in- 
volves the six partly filled 3d orbital on the Cu(II) sites 
and two lowest energy 4d orbitals on the Mo(IV) site. 
Besides the site energies of the orbitals and the transfer 
integrals between the Cu and Mo orbitals, we introduce 
on-site and extended range interactions between the vari- 
ous sites. On the Mo site, we consider all the three types 
of two electron integrals involving the Mo orbitals. Using 
this model, we compute the full eigenvalue spectrum, the 
transition dipoles with the ground state, spin and charge 
densities for all the states in each of the spin sectors by 
exactly solving the model in a spin conserving valence 
bond basis. Our studies show that the ground state is 
sixty four fold degenerate with an unpaired electron on 
each Cu site. This leads to the observed paramagnetic 
susceptibility prior to irradiation, of six non-interacting 
spin-1/2 objects per formula unit. The excited states to 
which there is nonzero transition dipole in each spin sec- 
tor are located at about 2.993 eV (for the values of the 
parameter we have chosen) in each of the spin sectors. 
The number of states at this energy as well as the transi- 
tion dipole moment for the transition arc highest for the 
S=3 states. There also exist very few states in the S=3 
sector below this energy, thereby introducing a bottle- 
neck for its internal conversion for non-radiative decay 
to the ground state. There may also occur changes in 
the geometry of the complex in the excited state which 
renders a low Frank-Condon factor for radiative decay. 
Thus, the observed post radiation magnetic behaviour of 
the complex with 75% of the molecules in the S=3 state 
can be rationalized by our model. 
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